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The aim of the present work is the investigation of solutions to the system of the Eulcr 
equations, including the systems with the right-hand sides describing various interior forces: 

dtp + div (pV) = 0, (1) 

d t (pV) + (pV, V) V + Vp = pf (x, t, V, p, S), (2) 

d t S + (V, VS) = (3) 

with the state equation p = e s p 7 , 7 = const > 1. We suppose that f is the smooth function of 
all its arguments. Here p(t, x), V(t, x), S(t, x) are the components of solution, corresponding 
to the density, the velocity and the entropy, given in K x K", n > 1. Equations (1-3) describe 
the balance of mass, momentum and entropy, correspondingly. 
Set the Cauchy problem for (1 - 3): 

p(0, x) = p (x) > 0, V(0, x) = V (x), S(Q, x) = 5 (x). (4) 

We deal with the classical solutions to system (1-3) with the density so quickly decreasing 
as |x| — > 00, to guarantee the convergence of the integral J p|x| 2 dx (so called solutions with 

R™ 

the finite momentum of inertia) . 

As well known, the solution to the Cauchy problem for system (1-3) may lose the initial 
smoothness in a finite time, sometimes there is a possibility to estimate the time of singularity 
formation from above (see, f.e., Q and references therein). Moreover, in the case more 
frequently investigated f = the compactly supported initial data are sufficient for the 
further singularity formation, (f.e., 0). 

At the same time it is interesting that there are some nontrivial classes of globally smooth 
solutions. Note that if f = 0, the components of such solutions do not belong to the Sobolev 
class. 

The paper is organized as follows. Firstly supposing the existence of globally in time 
smooth solution having concrete properties (denote it Uq), we show that if we choose some 
initial data (CDi) close in the Sobolev norm to initial data of such solution (CDo), then the 
corresponding solution to the Cauchy problem (E/i) occurs globally smooth as well. Then 
we show that the solutions Uq with such kind of properties exists, moreover, at n = 2 we 
construct some of them. 

The present work may be considered as a continuation of papers ||, [Q, and their 
generalization in some sense. 

Martially supported by RFBR grant 00-02-16337 
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1 Symmetrization and the result on a local-in-time smooth- 
ness 

The result on a local in time existence of smooth solution U(i, x) for symmetric hyperbolic 
systems, i.e. systems of the form 

f)J j " f)\ j 

A°(t, x, U)— + ]T A k (t, x, U)^ = 5 (i, x, U), (1.1) 
fc=i 

where the matrices A J (t, x, U) are symmetric, and additionally the matrix A°(t,x, U) is 
positive definite, is well known the particular cases in [Q, ||, ||)). 

Namely, if the matrices A 3 (t, x, p), g(t, x, p) depending smoothly on their arguments, 
have continuous and bounded derivatives with respect to the variables (x, p) up to order 
rn+1 under bounded p, the initial data Uo = U(0, x) and the function g(t, x, 0) belong to the 
Sobolev class H m (W l ), m > l + n/2, for any fixed t > 0, then locally in time the correspond- 
ing Cauchy problem has a unique solution from the class n] =o C J ([0, T); B m ~ 3 (R n )), T > 0. 
Moreover, 

Jm sup(||U||i- + ||V x U||i») = +oo, (1.2) 
T = T(U ) and lim T(U ) = +oo. 

||Uo|| H m-Q 

(Note that less rigid requirements may be imposed on the initial data for the existence 
of smooth solution locally in time, f.e. local (in space), but uniform belonging to H m (M. n ). 
1,0). 

In many problems having the origin in physics the coefficients of system (1.1) depend 
only on solution, this simplifies significantly the formulation of the result. In the case for 
the local in time existence of the Cauchy problem in the Sobolev class described above, 
it is sufficient, beside of the coefficients smoothness, to require the implementation of the 
condition g(0) = 0. 

As we deals with the solutions to system (1-3) such that infp = mfp = 0, we have to 
use the symmetrizartion firstly proposed in pi . 

Involving the variable II = k (pi 2) 2 ~< , k — — - — , we get the symmetric form as follows: 

7-1 

exp(-)(dt + V, V)n + ^— exp(-)nd*0 V = 0, (1.3) 

7 2 7 

(ft+(V,V))V+^=^ex P (-)nvn = f 1 (t,x,n,V,5) = f(i J x,e-T ,V,S), (1.4) 

2 7 \nj 

(d t + (V,V))S = 0. (1.5) 

Denote U = (II, V, S). As an immediate corollary of the general result on the symmetric 
hyperbolic systems we obtain the following 

Theorem 1. 1 Let the initial data Uo = (IIo, Vo, So), belong to the class H m (M. n ), m > 
1 + n/2. Suppose the function fi(t, x, U)) have continuous and bounded derivatives with 
respect to space variables and the solution components up to order m+1 for bounded U, and 
fl(t, x, 0) € H m (W l ) for any fixed t > 0. 
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Then the Cauchy problem for system (1-3) has locally in time a unique solution such 
that 

(II, V,S,)G n} =o Ci([0,T);H m -i(R n ), 

moreover 

lim T = +00. 

||U || H m-0 

2 The interior solution 

We call an interior solution globally in time smooth solution (p(i, x), V(t, x), S(t, x)) to sys- 
tem (1—3) having the following property: another solution (p, V, S) to the Cauchy problem 
(4) having sufficiently small norm 

||04 7 - 1)/2 (x) - ^(o.x), V (x) - V(0,x), 5 (x) - S(0,x)\\ Hm{Mn) 
is smooth globally in time as well, moreover 

(p(7-D/2 _ j5 (7-i)/2 ) v _ v, s- S) G r^ =0 C s ([0,oo);H m - i (R n )), m > l + n/2. 
Note that the trivial solution is not interior if f = 0. 

The set of interior solutions is not empty. In the paper |Q (the generalization of ||) for 
f = it is shown that the solution to system (1-3) (0, V(t, x), const) is interior, if V(t, x) 
is the solution to equation 9 t V + (V, V) V = such that SpecDV(0, x) is separated from 
the semi-axis K_ and DV(0,x) e L°°{R n ), D 2 V(0,x) e if m_1 (R™). (We denote D k the 
vector of all spatial derivatives of order k.) 

In particular, as it is noted in Q| in the role of such solution may be played by the 
solution with linear profile of velocity V = A(t)r, where r is the radius-vector of point, A(t) 
is the matrix with the coefficients depending on time such that Spec^4(0) ^ R_. The solution 
may be found explicitly. 

Note that the given interior solution (and the close solution as well) does not belong to 
the Sobolev space, since the velocity components grow at infinity. The result is clear from 
the physical point of view: the velocity field having the positive divergency " spreads out" 
the initially concentrated sufficiently small mass, that prevent the singularity formation. 

But the given solution with the zero density is not physical. It follows, in particular, in 
the case of the velocity Jacobian having initially negative spectrum at some value of x , the 
appearance during a finite time of points with infinite negative divergency. If the solution 
with non-zero velocity, initially close to the described above, would behave like this, it would 
follow (according to (1)) the infinite value of density in some points. 

The question arises: can one construct the interior solution with the density not close 
to zero in the Sobolev norm? Or to be interior the solution must describe the scattering of 
sufficiently small mass? 

If from some reasons we succeed to find the velocity field, then from the linear with 
respect to p and S equations (1), (3) we can find the solution to the Cauchy problem (4), 
moreover, if the velocity field occurs smooth, then the rest of the components are smooth. 

It needs to note, that the values p and p can become unbounded in some point of 
trajectory in a finite time T if and only if J Q T divVdt = —00 on the trajectory (see (1)). 

Before formulating the theorem we do some transformation under supposition that 
(B, V = A(i)r, S) is a globally in time smooth solution to system (1.3-1.5). 
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According to (1.3-1.5) the vector-function u = U— U (here u = (tt, v, s), U = (II, V, 5), U = 
(fl, V, S)) satisfies the system of equations 

(fit + V, V)tt + (v, V)tt + (tt + n)divV = 

-(v,V)n-^7rdivV, (2.1) 

{d t + (V, V))v + (v, V)v + ^— - exp(-)(7r + n)V7r = 

2 7 

-(v, v)v - ^— i exp(-)(7r + n)vn + exp(-)nvn+ (2.2) 

27 27 

+fi (t, x, (V + v), (S + s) , (S + tt)) - fi (t, x, v, s, n) , 

(d t + V,V)s+(v,V)s = -(v,V)S. (2.3) 

Further, following to we carry out the nondegenerate change of variables such that 
the infinite semi-axis of time turns out to semi-interval [0, Coo), to obtain the symmetric 
hyperbolic system with the coefficients allowing to apply the theorem on a local existence 
of smooth solution on the interval [0, cr*), tr* < <7oo. 

So, let (ti,y) :— (t,A(t)x) be the new variables (A(t) = \\a,ij(t)\\ is a quadratic (n x n) 
matrix). Then V x = A*\7 y: div x V = div y ^4V, d tl — d t + ArX7 y . Choose the nonnegative 
function A(f) such that the integral J Q A(r)dr = < oo converges and set <r(t\) — 

Jq 1 A(r)dr. In that way, the semi-infinite axis of time goes to the semi-interval [0, Coo). 
Further, involve the variables 

W = A(t)vA _1 (t), P = A- 9 (£)tt, 

the constant q will be defined below. 

After the transformation we get the system 

{d c + (W, V y ))P + -^(P + P)div y W = 

-(W,V y )P-PQ 1 (a), (2.4) 

(d a + (W, V„))W + + P)V,P = 

_2_Ivi;(5, <r)(P + P)V,P + ^)PV,P - WQ 2 (<r) + G, (2.5) 

(S ff + W,V v )a = -(W,V v )5. (2.6) 

Here we denote 

Qi(t(v)) = A(f)" 1 (^tr A(t) + <?(mA(i))'), 
Q 2 {t{a)) = \(t)- l ({\n\{t))'E + AiWA- 1 )'® + E)), 

9(S,a)=aq?(-)R(a)B(p), 

7 
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B(t(a)) = A{t)A*{t){detA{t))- 2 ' n 

R(t{a)) = A 2 «- 2 (t)(dctA(t)) 2 /", 

G = A- 2 (t)A(t)(f!(i, A(t)x, \- q (t)(P + P), A(t) J 4 _1 (t)(W + W), (S + «))- 

fi(t, A(i)x, A- g (*)P, A^JA-^tJW, 5)), 

taking into account that t — t(a). 

Note that B is a bounded invcrtiblc matrix. 

After multiplying (2.5) by ^~ 1 (S, a) the system become symmetric hyperbolic, however, 
generally speaking, ^^ 1 (S, a) — > as t — > oo or will be unbounded owing to the behaviour 
of -R(o-O)). 

In another words, we can increase the time of existence of smooth solution to system 
(1-3) as long as we wish due to choosing the initial data small in the Sobolcv H m norm, 
but we cannot guarantee the existence of the solution during the infinite time. Let 

D "-{w)"-\wT ■?<*--- 

Denote Cj~(M. n ) the space of continuous functions having continuous bounded in 1" 
derivatives up to order r = 0,1, .... Further, denote fv the matrix 

Now formulate the theorem taking into account all denotation involved. 

S 2 

Theorem 2. 1 Let the function fi(f, x, II, V , S) =f(t,x,e~~ >V, £) have the deriva- 

tives with respect to all components of solution up to order m + 1, continuous and bounded 
under bounded (II, V, S), moreover f\(t, x, 0, 0, 0) G H m (R n ) at any fixed t > 0. Suppose 
system (1-3) has the globally in time smooth solution U = (p,V,S) with the linear profile 
of velocity V = A(t)r such that 

a) p^(0,x) G n^+^OR™); DS(0,x) G n™ ^'(K"); 

b) = dct A(t) >0 fort>t >Q; 

and there exist a smooth real-valued function X(t), a constant q and a skew symmetric matrix 
with the real-valued coefficients U<f,(t) with the following properties: 

c) / t + °° \{r)dr < oo, 

rf;it oo A«(T)^/"( T)dT<0O; 

e) functions 

Qi(t)R(t), (2.7) 
Q 2 (t) - X-^AtyfyA- 1 ® - U^t)), (2.8) 
(Infl(i))' (2.9) 

are bounded on M + x R n . 
Then the solution U is interior. 
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Proof. At first we consider the case of R(a) such that < r\ = const < R(<r) < r 2 = 
const < oo. If the other conditions of Theorem 2.1 hold, then we can prove Theorem 1.1 
immediately. The case is realized, for example, at 7 = — + 1, f = 0, in the situation, 

considered in § {trA(t) ~ f, £(t) ~ t~ n (t -> 00), A(t) = (1 + t)" 2 , g = 



In the general case the proof is more complicated. 

Suppose a £ [0, cr») and (X, Y, Z) : i?" — * i? x i? n x is a vector-function from 
Following ||, involve the norm 

[A,Y,Z] 2 (<t) := /" (X 2 + Z 2 +Y**- 1 (5, ( r)Y) ( iy. 

For p £ N define 

e p(°) = \H { D ° p > ° aw > ^V), 

\ct\=p 

m 
p=0 

If we should succeed to show that F m (a) is bounded on [0,<7oo] function, the theorem 
would be proved, as the impossibility to prolong the solution is connected with the going to 
the infinity of L^— norms of the solution itself or its gradient (see (1.2)). 

Showing that H m — norm of the solution admits a majorization by means of some func- 
tion, bounded on [0, <7oo], we shall show that cr* = tJoo and thus the theorem will be proved. 

After involving of the new variable II, which is necessary for the symmetriization of the 
system (see above), globally smooth solution U corresponds to the globally smooth solution 
(n, V, S) with the initial data (flo, Vo, So)- After the changing of variables it corresponds 
to the solution (P, W, S) with initial data (Pq, W , Sq). The last solution is a smooth 
solution to system (2.4-2.6) on [0, o"oo], and therefore its H m — norm is bounded on the 
segment. 



Compute 

^ = E [/ {D a PD a 8 a Pdy- 



D a sD a d a sdy+ 



J D a W^-\S, a)D a d a Wdy+ 

Rn 

[ (D a W)*( d ^ - -^- 1 (S,a)d a S)D a Wdy] 

J da 7 
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h+h+h+ h 

(further it is supposed the summation over repeated indices, the integrals Ik, k = 1,2,3,4, 
are numbered in consecutive order). 

Estimate all integrals. Below we denote by a the positive constants, depending only 
on initial data. Begin from I4. Before all we stress that according to (2.9) R'(a)/R(a) is 
bounded as a — > a 00l therefore 

f^lf'h = 1 - '-'WW*"* + oct-w.))! < 

Thus, 

h < c 2 (l + HWVyslloo + ||WV y 5||oo)-B P < (C2 + c 3 P 1/2 P m + aR^F^Ep. 
Further, 

h=- J D a PD a ((W, V y )P + + P)div y W + (W, V y )P + PQ^dy. 

Rfl 

The integral from the member of higher order 

7-1 
2 

been integrated by parts can be reduced to 



D a PD a ((W, V y )P + \- Pdiv y W) 



- J I i|D Q P|div y W + D a P ((W, V y (P Q P)) - D a (W, V y P))| dy- 
- y P a P{(V y P,P"W) + (Pdiv y P Q W-P Q (Pdiv y W))}dy :=J n . 
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2 

Under sign of the integral there is the sum of members of the form d a Uid l U2d a+1 ~ l U 2 , < 
I < a (by Uj, j — 1,2, we mean the corresponding components of the solution. At first let 
a ^ 0, 1. According to the Galiardo-Nirinberg inequality we have 

WU 3 \^<C\DU,^^\D a U^ 

for pi = 2^5jk For I ^ 1, a it is true that — + - — = \. From the Holder inequality it 
follows that 

J \d a U 1 d l U 2 d a+1 - l U 2 \dy < 

Rn 

||p a ^i||2||a't/ 2 9 a+1 -'[/ 2 || 2 < 
\\D a u 1 \\ 2 \\d l u 2 \\ m \\d a+1 - l u 2 \\ Pa _ l+1 <c\\D a u 1 \\ 2 \\D a u 2 h\\D a u 2 \\ <x . 

For the other values of a and I the last inequality is evident. 
Thus, 

hi < cHDPPlladlDPPllallVyWHoo + ||P P W|| 2 || V y P|U) < cR^F^Ep. 
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(one can find the details of the inequalities applying in the book 
Further, according to condition (2.7) of the Theorem 



D a PD a (WVy-P)dy 



< c 5 R 1 / 2 F m , 



J D a PD a (PQ 1 )dy 



< ceQiRF m < cjF„ 



The integral 1 2 can be estimated analogously: 

\h\ < c 8 R^ 2 F^ 2 E q + CgR^Fr, 

At last, 



h = J ( J D Q )*W*- 1 (5,CT){(W,V y ) J D Q W- J D Q ((W,V y W)}dy+ 



(D a W)*V-' L (S, cr)L» a W)div y W - WV y S (ywjT^S, cr)L> Q W j dy+ 



7-1 



{(D a W)*^- 1 (S, a)D a (PV(S, <r)V y P + PV(S, <r)V y P- 

(P + P)V{S, <j)V y P + PV(S, a)V y P) } dy+ 
J (D a W)** -1 (S, (WQ 2 + G) dy. 

The first and second integrals can be estimated from above by values 

cioir^wiiliivwiioo + ciiir^^ 

and the third one by the value 

ci3\\D p WU\\DPp\\ 2 + \\DP.s\\ 2 + \\DVP\\ 2 + ||^- 1 V y 5|| 2 )(l + \\P\U + ||V y P|| 

||V y fl + Hoc H^Hoo + l|V y P||oc + \\VyS\\ 00 ) P+1 < 

c 13 R^( Cli + F^/ 2 )E P + c 15 R^ 2 (c 16 + F^I^El'*, 
and, at last, the fourth one taking into account conditions (2.8) by the value 

cnE p . 

In the last case note that 

G = A- 1 (t)A(t)f v (t,x,n,V + flv,5)>4- 1 (t), Oe (0,1), 



8 



J(D a W)*^- 1 (S, a)D a (U c k(t)W)dy = 0, 

Rr, 

as (w, U^(t)w) = for any vector w G R n . 
Uniting all the estimates we get 

E' p < c 18 E p + c 19 R 1/2 E p + c 20 R 1/2 F^ 2 E p + c 21 R 1/2 F m E p + 
c 22 R 1,2 F^' 2 E p + c^R^F^^El' 2 , 

whence it follows after summation over p that 

F' m < c 2A F m + c 2b R^ 2 {FU 2 +F m + FU 2 + F^' 2 + F^ 2 ). 
Set A m (cr) = cxp(— C24<r)F m . Then 

< c 26 (AV 2 + A m + K% 2 + A^+ 2)/2 + A^ /2 )R 1/2 , (2.10) 



where the constant c 2 q depends on c 2 4, c 2 §, a a 
We set 

g 



5>o 



dg 



g 1 / 2 + g + g 3 / 2 + g('P+ 2 )/ 2 + g(p+ 3 )/ 2 ' 



The integral diverges in the zero, therefore 0(0) = — oo. 
Integrating inequalities (2.10) over a we obtain that 

a 

0(A m (a)) <6(A m (< )) + c 26 J R^^da, 

ff(to) 

moreover, as follows from (2.8), the integral in the right-hand side of the last inequality 
converges as a — > (Too to the constant C, depending only on the initial data. 

Choosing A m (to), (and together with the value the H m — norm of the initial data) suf- 
ficiently small, one can obtain that 0(A m (io)) + C is later then 0(+oo), it signifies that 
A m (a) and F m {a) are bounded from above for all a S [0,Uoo) and er» = a^. So, Theorem 
2.1 is proved. 



Remark. The solutions with the linear profile of velocity, satisfying the theorem condi- 
tions, exist. In the case described in (| / = 0, 7 < 1+f , A(t) = (J5+A(0)t) -1 A(0), SpecA(O) i 
R_, trA(t) - f,detA(i) ~ i, Q x = 0,Q 2 ~2 (i((trA _1 )£J - A -1 ) , as i -> 00. Here 

\=(l + t)- 2 ,q=^, 17^=0. 



We shall get two corollaries from Theorem 2.1, basing on which we can assert that the 
solutions which will be constructed in the next section are interiors. 

Corollary 2. 1 Let f = 0. If system (1-3) has the globally smooth in time solution with the 
linear profile of velocity, described in the statement of Theorem 2.1 and A(t) ~ jE, t — » 00, 
where d is a positive constant, then the solution is interior. 
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Proof. Choose X(t) = t s+1 , q > | — 2 (s+\) , U$(t) = 0. In that way all conditions of 
Theorem 2.1 are satisfied. 

Corollary 2. 2 Lei f = LV, where L is a matrix with the smooth coefficients such that 
AqLAq 1 — U\{t) = —(J>E, where Aq = 5E + U% is a matrix with the positive determinant, the 
matrices U\ and U% are skew symmetric, 5 is a constant, fi is a positive constant. If system 
(1-3) has the globally smooth in time solution with the linear profile of velocity described in 
the statement of Theorem 2.1 and A(t) ~ jAq, t ~ * °Oj then the solution is interior. 

Proof. To verify condition (2.8) one has to show the boundedness of the value 

A -1 (t)((ln X(t))'E + A(t)E - A(t)LA -1 (t) - A-\t)A\t) - U+) ~ 

\- 1 (t)((ln\(t))'E+ -((8+ l)E + U 2 ) - Ui(jb)+nE-U4,),t-KX>. 

Choose X(t) — t s+1 exp{— /it}. Condition (2.7) is satisfied at q > §, condition (2.9) can be 
verified elementary. It needs to choose the matrix U^t) = j\Ji — Ui(t). 



3 Constructing of interior solutions 

Limit ourselves to the important case of n = 2 and set f = Lv, where the matrix L = 

( j 1 ^ J , fi is a nonnegative constant, I is an arbitrary constant. Thus, we in the 

simplest way describe the Coriolis force and the Rayleigh friction, for example, in the me- 
teorological model under neglecting the vertical processes and the Earth curvature. 

It is easy to verify that on the smooth solutions to system (1-3), so quickly vanish- 
ing at infinity to assure the convergency of all integrals involved, there are the follow- 
ing conservation quantities: the mass m = J pdx = const, for /i — the total en- 

R 2 

ergy E = j ( p ^ ^ H — 2_) dx = E^(t) + E p (t) — const and the momentum of inertia 

R 2 7 ' 

J = J /o((Vj_,r)+ 5l x | 2 )^ x = const. The conservation laws are true also for the solutions 

R 2 

to system (1-3) with the density (and the pressure) quickly vanishing as |x| — > 00, and the 
velocity components may even increase. 

We consider another integral functionals, characterizing the average properties of solu- 
tions, namely 

G(t) = ±f p\r\ 2 dx, Fi(t) = JiV^Apdx, 1 = 1,2, 

R 2 R 2 

where X x = r = (x, y), X 2 = r± = (y, -x). 
Note that G(t) > 0. 

On smooth solutions to system (1 - 3) the following relations take place |H| 

G'(t)=F 1 (t), (3.1) 
F^t)=lF 1 (t)- f xF 2 (t), (3.2) 
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F{(t) = 2( 7 - l)E p (t) + 2E k (t) - lF 2 (t) - ftF^t), (3.3) 
E'(t) = -2»E k (t). 

Remark once more that at (i = the quantity E(t) of total energy is constant. The functions 
E p (t) and E k (t) generally cannot be expressed through G(t), F\(t), F2(t). 

We try to choose the form of velocity field in such way that from system (3.1 - 3.3) one 
could find it directly. Suppose, for example, 

V = A(t)r, (3.4) 

where A(t) is a (2 2) matrix , whose time-depending coefficients we have to find. 
1) We get the simplest result if we choose the matrix A(t) such that 

V = a(t)r + /3(t)r ±) 
It is easy to sec that in the case we have 

= 2a(t)G(t), F 2 (t) = 20(t)G(t), 

E k {t) = {a 2 {t)+f3 2 {t))G{t). 
It is not difficult to obtain from the equation 

d t p+ (V.Vp) + jpdivV = 0, (3.5) 

which is a corollary of system (1-3) and the state equation, that 

Ep (t)=E p (0)G^-\0)-^ T ^. 

In that way, all functions involved in system (3.1 - 3.3) are expressed through G(t), a(t), f3(t). 
For the sake of convenience we denote Gi(t) = 1/G{t) and get the system of equations 

Gi(t) = -2a(t)Gi(i), (3.5) 

(3'(t) = a(t)(l-2(3(t))-^(t), (3-6) 
a'(t) = -a 2 (t) + (3 2 (t) - i/3(t) - m {t) + ( 7 - l)E p (0)G{-~< (O)Gj(t). (3.7) 

The solutions of the system of ordinary differential equations are smooth if they are can 
be prolonged up to all axis of time. The non-prolongation is connected with the escape of 
the components of solution to infinity during a finite time. 

Note that for 7 > 1 the quantities a(t) and (3(t) are bounded, this follows immediately 
from the expression for the total energy of the system: 

E(t) = (a(t) 2 +(3{tf)G{t) +E p (0)G^\0)-^ r) < E(0), 

therefore 

a 2 (t) + )3 2 (t) < £7(0)Gi(t) - E p (0)G\-~< {0)G~[ (t) < +00, 

thus the quantities of density and pressure remain bounded for all solution we want to 
construct. 
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At [i = system (3.5-3.7) can be integrated explicitly: 



a(Gi) = ±<JkGJ - C 2 G\ + {E- /C)Gi - I 2 /4, 

m=CG 1 (t)+l/2, (3.8) 

f Gl(t) _dGj_ 
J Gl{0) 2G ia (Gi) 

with the constants G = # = (a 2 (0) + G 2 G?(0) — (E — ZG)Gi(O) + / 2 /4)/G*7(0). 

Due to (3.8) one can reduce the number of equations in the system and consider it on 
the phase plane (Gi, a), G\ > 0. There exists the unique singular point: at I =/= it is a 
center situated on the axis a — 0, at I ~ it is a complex equilibrium in the origin, the 
trajectories form the elliptic saddle point. In the last case the time of movement from any 
point of the phase plane to the origin is infinite. The equilibria, at I ^ lying in the plane 
G\ = , namely, the points a = [3 = and a = 0, (3 = I, are the centers under consideration 
in the plane (a, (3). 

At [i 7^ one needs to integrate the system (3.5-3.7) numerically However, the conclu- 
sion on the behaviour of the solution components at infinity we can do analytically as well. 
At /i > the system has a unique stable equilibrium in the origin (at I = it is a knot, at 
I ^ it is a focus.) Therefore we can find the solution as a formal asymptotic series by the 
negative powers of t. Limit ourselves to writing out the leading terms: 

a(t) ~ M p , (3(t) ~ M r , Gi(t) ~ cii 9 , a 1 ,b 1 ,c 1 = const ^ 0. 

From (3.5) we have immediately that qc\t q ~ x — —2a\C\t v + q , therefore?? = —1, q = —2a\. 

From (3.6) we obtain that rb x t r ' Y ~ la x t v - MM r - 2a 1 b 1 tP +r . If I ^ 0, fi ^ 0, then 
r = p = —1, moreover, la\ — \xb\. 

From (3.7) we get analogously that pa^ 1 a\t 2 P + b\t 2r - lb x t r - /j,ait p + Kc\V q , 

where it is denoted for short K = (7 - l) J B p (0)G^ 7 (0), that is -aii~ 2 aft" 2 + 6 2 ^ 2 - 

Iht' 1 - fiaxt- 1 + Kc[P q . 

Let 7(7 < —1, that is the last term does not contain the senior degree. Then the rest of 
members of senior order must be eliminated, and — Z61 = fiai = fi 2 bi/l, therefore I = = 
in spite of the supposition. Consequently, q — — 1/7. 

So, at fi ^ 0, I ^ 

a(t) . _U-, m „ i t -., Gl(t) „ (<il±£f , _ „. 

At I = 0, due to the existence of the integral /3(i) = GGi(i) exp{— /it}, system (3.5 - 3.7) 
can be reduced to the nonautonomous system of two equations: 

Gi(*) = -2a(t)Gi(t), 

a'(t) = -a 2 (t) - /j,a(t) + C 2 cxp{-2 M t}G 2 (i) + KG'((t). (3.9) 

From (3.9) we have that pa^ 1 a\t 2 P - iicut? + C 2 c\t 2q exp{-2^} + Kc\V q , where, 

remember, p = —1, q = —2a\. To compensate the senior term —[lait^ 1 it is necessary that 
the equality q = — 1/7 holds. Hence, we find the coefficient c\. 
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Thus, at fj, ^ 0, I = 

a(t) ~ ^t- 1 , /?(*) - G (^) V7 t"V7 exp^i}, Gl - (^) ^ t" 1 /?, t _ oo . 
At p = I — 0, as follows from the explicit formulas, 

a(t) ~ (yG^jOj/E + tj 1 , ~ a 2 (t), = GGi(t), t -» oo. 

Thus, as follows from Corollaries 2.1 and 2.2, in the cases it = I = and /x > we 
have constructed the velocity field for the interior solution. More precisely, at p, = I = 
A(t) <~ ji?, i — > oo, Corollary 2.1 can be applied with 5 = 1. At > one can apply 
Corollary 2.2. In the case /x > 0, I — we have U\ = Ui = 0, S = If p, > 0, Z ^ 0, then 

As soon as a(t) and (3{t) are found, the rest of the components can be found elementary 

t t t 

p(t,\r\,((>) =exp(-2 J a(T)dT)p (\r\cxp(- J a(T)d,T),<j) + J /3(r)dr), 



t t 
S(t, |rU) = S (|r| exp(- y a(r)dr),0 + | /?(r)dr). 

o o 
t i t 

=exp(-27 J a(T)d,T)p (\r\exp(- J a(T)d,T),(f> + J /3(r)dr). 



From (2) and (3.7) we get that on the classical solution of the initial system must be satisfied 
the relation 

Vp = -( 7 -l)Gl- 7 (0)£? p (0)G7(t)pr. 

Hence it follows that the components of the initial data po po must be axisymmetric and 
compatible, i.e. connected as follows: 

Vpo = -(7-l)Gi(0)E p (0)por. 

However, in spite of the components of the pressure and the density must vanish as 
|x| — > oo, the entropy must even increase (remember that the conditions to Theorem 2.1 
require only the boundedness of the entropy gradient). However, one can choose 

1 2a 1 

^° = /-, . i i?w a = const > 3, po 



(l + |r| 2 )"' ' ^ (7- l)Gi (0)^(0) (l + |r|2)o+i' 

then S = const + (0(7 - 1) + 7) ln(l + |r| 2 ). 

Remark. In the case G(0) ^ the density cannot be compactly supported (in contrast 
G(0) = 0). Really, since p — tt^- 1 , p = ir-t- 1 e 1 , then due to the compatibility conditions 
7rV7r ~ const ■ e~~, |x| — > c — 0, where c is a point of the support of n. Therefore, for C 1 — 
smooth 7r it occurs that S — > +00, |x| — > c— 0, and we cannot choose any smooth initial data. 
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It is interesting that if one requires only C°— smoothness of it (tt <~ const- (c— Ixl) 1 / 2 , |x| — » 
c — 0, 7T = 0, |x| > c > —0), the condition may be fulfilled. Moreover, p and p will be of the 
C 1 class of smoothness, however, neither the theorem on the local in time existence of the 
smooth solution, no Theorem 2.1 can be applied. 



2) To consider the velocity field (3.4) with the matrix 

AW- 



ait) b{t) 
c(t) d(t) 



involve the functions 



G x{t) = ^ J px 2 dx, G y {t) = i J py 2 dx, G xy (t) = i j pxydx. 

R 2 R 2 R 2 

Then involve the auxiliary variables 
Gi(t) = G x {t)^ [l+1)/2 {t), G 2 (t) = G y (i)A-^ +1 )/ 2 (i), G 3 (t) = G xy A-(T +1 )/ 2 (i), (3.10) 

where A(t) = G x (t)G v it) — G xy (t) is a positive function on solutions to system (1-3). Note 
that the behaviour of A(i) is governed by the equation 

A'(t) = 2(o(t) + d(t))A(t), (3.11) 

and a potential energy -E p (i) is connected with A(t) as follows: 

E p (t) = ^ p (0)A^- 1 )/ 2 (0)A(-T+ 1 )/ 2 (f). 

To find Gi, G2, G3, and also the elements of the matrix Ait) we get the system of 
equations 

G[(t) = ((1 - 7 )a(t) - (1 + 7)d(*))Gi(t) + 2b(t)G 3 (t), 
G£(t) = ((1 - 7)d(t) - (1 + l)ait))G 2 (t) + 2c(*)G 3 (i), 
G 3 (i) = c(t)Gi(t) + 6(i)G 2 (t) - 7(a(*) + d(t))G 3 (t), 
a'(t) = -a 2 (t) - &(i)c(i) + ic(t) - /xa(t) + i^iG 2 (t), (3.12) 
6'(t) = -6(t)(a(t) + d(i)) + Jd(t) - /i6(i) - #iG 3 (i), 
c'it) = -c(i)(a(i) + d(t)) - io(i) - //c(t) - #iG 3 (t), 
d'(i) = -d 2 (t) - b(t)c(t) + lb(t) - fid(t) + Kidit), 
with K x = ^(^(OjAfr- 1 )/ 2 ^)). 

Proposition 3. 1 In i/ie case p = I = any solution to system (1-3) with a linear profile 
of velocity, satisfying the condition ) of Theorem 2.1, is interior. 



14 



Proof. Show that at [i = I = there is the asymptotics A(t) ~ jE, t — » oo, that 
is according to Corollary 2.1 the corresponding solution with a linear profile of velocity is 
interior. 

Go to the new variables a^t) = a(t) - d(t),bi(t) = b(t) + c(t),a(t) = b(t) - c(t), di{t) = 
a(t) + d{t),G 4 {t) = G 1 (t) + G 2 (t),G 5 (t) = G 1 (t)-G 2 (t). In variables a u b u c u d u G 3 , G 4 , G 5 
system (3.12) has a form: 





a[(t) = -ai(t)di(t) - KiGsit), 


(3.13) 




b' 1 (t) = -b 1 (t)d 1 (t)-2K 1 G 3 (t), 


(3.14) 




c[(t) = -ci(t)di(t), 


(3.15) 


d' 1 (t) = - 


-i(a?(t) + d?(t)) - |(6?(t) - c?(i)) + JfiG 4 (t), 


(3.16) 


G'S) - 


= - 7 di(i)G 3 (i) + ^i(t)G 4 (t) - ^ci(t)G B (t), 


(3.17) 




= - 7 di(t)G 4 (t) + ai(t)G B (t) + 26i(t)G 3 (t), 


(3.18) 


G'S) 


= - 7 di(t)G B (t) + ai(t))G 4 (t) - 2ci(t)G 3 (t). 


(3.19) 



Let as t — > oo the asymptotics of functions involved in the system be the following: 
oi(t) - Lit' 1 , 6i(t) - L 2 t h , Cl (t) ~ i 3 t' 3 ,di(t) - L±t l \G-i{t) ~ iVi', G 4 (t) - M^ 1 , 
Gs(i) ~ M 2 t P2 , where Lj, i = 1,2, 3, Mj, j = 1,2, iV are some constants not equal to zero. 
Note that p 2 < p\ and in virtue of A > the estimate q < pi holds. 

From (3.15) we get immediately that 

Z 3 L 3 *' 3-1 = -L 3 L 4 t h+h , 

hence U = — 1, L 4 = — i 3 . 

From (3.13) we have taking it into account 

hL^ 1 - 1 = -L^t 11 - 1 - 2K 1 M 2 t P2 . 

The following variants are possible: 

Pi < h — 1, hence Zi = — L 4 , i.e. l\ = I3, 
p 2 = h — 1, hence 

_,,_*-«!. (3.20) 

From (3.14) we have analogously 

hL^- 1 = -L 2 L^- X - 2K\Nt q . 

There are the variants: 

q < l 2 — 1, hence l 2 = —L 4 , i.e. l 2 = l 3 , 
q = l 2 — 1, hence 

, 2 = _ l4 _^ =i3 _^. (3. 21) 

From (3.16) we get 



-Lit- 2 = ~{Llt 2h + L 2 t- 2 ) - ^(L 2 t 2h - L 2 t 2h ) + K 1 Mit Pl . (3.22) 
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If k < — 1, i = 1,2, 3, then at t —* oo 



a(t) ~ + L4*- 1 ), 6(0 ~ ^2*' 2 + L 3 t l % 

c(t) ~ i(i 2 i i2 - L 3 t ,s ), d(t) ~ ^t -1 - 

that is - ^E, L 4 = const > 0. 

We shall show below that there are not others situations. 

From (3.22) it follows the impossibility of the situation l 3 = —l,l\ < — l,l 2 < — l,Pi < 
-2. Really, in the case l 3 = -L 4 = -1-y/l + L\ < -1 or l 3 = -L A = -1- y/1 + L\ + 2K 1 M 1 < 
— 1, it contradicts to the supposition. 

Now let l 3 > —1. Then from (3.22) l\ — l 3 or (and) l 2 — h- Suppose, for example, that 
li = l 3 . Then p 2 < l\ — 1, from (3.19) it follows that 

p 2 M 2 t P2 ~ x = -jUM^ 2 - 1 + L 1 M 1 t Pl+h + 2L 2 Nt q+l \ (3.23) 

and therefore pi + h = q + l 2 , h < h, i.e. l 2 > — 1. Besides, 

LiMi = -2L 2 N. (3.24) 

From (3.17) we have 

qNt"- 1 = -jLiNt"- 1 + ^L 2 M 1 t Pl+l2 - ^L 3 M 2 t P2+h . (3.25) 

As soon as q — 1 < pi — 1, and p\ — l < pi + l 2 , then pi+l 2 = p 2 + 1 3 < p\ + 1 3 , and therefore 
h < '3, and as Zi = Z 3 , then l\=l 2 = l 3 and pi = q. Besides, 

L 2 M X = L 3 M 2 . (3.26) 

Further, from (3.18) we obtain 

Vl M x t Vl - x = -jL 4 M 1 t Pl - 1 + L 1 M 2 t P2+h - 2L 3 Nt q+h . (3.27) 
As q + l 3 > pi - 1, then p 2 + h = q + l 3 , p 2 = q, 

L 2 M X = 2L 3 N. (3.28) 

From (3.24) and (3.28) we have f| = -j*, and from (3.26) we obtain jfe = , i.e. 
L| = — L3, L2 = L 3 = in spite of the supposition. 

Suppose now that l 2 = l 3 . Then q < l 2 — 1, p 2 < li — 1. From (3.17) we have that 
1 < Z2 + Pi, since Z 2 > — 1, and therefore l 2 + pi = l 3 +p 2 , p\ = p 2 , L 2 M\ = L 3 M 2 . 

Taking this into account from (3.23) we have that pi + h = Q + h < Pi + h, h < 
l 2 , L 2 M X = -2L 2 N. 

From (3.27) we get p 2 + h = q + h < Pi + h, P2 = q, h = h, L\M 2 = -2L 3 N. In that 
way, we obtain the contradiction analogous to the previous one. 

Now let l 3 < -1, h > -1 and (or) l 2 > -1. Then from (3.22) it follows that Pl = 2h 
and (or) pi — 2l 2 (p 1 > —2). For example, if pi — 2/ l7 then from (3.24) we get p 2 — 1 < 
Pi+h, Pi+h = q + h < Pi+h, h < h- From (3.25)p 2 + ' 3 < h+Pi > h+Pi > 
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therefore q — 1 = fa + Pi- But q—1 < pi — 1, therefore fa < — 1, and Z 2 < — 1 in spite of the 
supposition. 

If pi = 2fa, then from (3.25) we get I3+P2 < Pi — 1, therefore pi + Z2 >= P2 + £3, <? — 1 = 
pi + I2, fa < fa. But q — 1 < pi — 1, therefore Z2 < — 1 in spite of the supposition. 

So, it remains the unique possibility: fa < — 1, Zi < —1 and (or) Z2 < — 1- In the 
follows from (3.22) 

Z 3 - -1 ± V 1 - (^1 + ^2 - W - 2^iMi), (3.29) 

where <5j = 1, if Zj = —1, and 5 j = otherwise, i = 1, 2, 3, (54 = 1, if p 2 = —2 and £ 4 = 
otherwise. Hence 

5iL\ + S 2 L 2 2 - 5 3 L 2 3 - 28±K X M X < 1, (3.30) 

if inequality (3.30) is strict, then fa < —1. 

We consider this case, that is fa < — 1, fa = — 1 and (or) fa = — 1. 

If Z x = -1 ^ Z 3 , then p 2 = fa - 1 = -2. If Z 2 = -1 7^ Z 3 , then g = fa - 1 = -2. In that 
way, in any of these cases we have pi = —2. 

At fa = -1, Z 2 < -1, < -1, from (3.23) we obtain p 2 = -7X4 + 

^ --<> + ■*). (3.31) 

from (3.27) we get Pl = -7L4 + 

L1M2 
Mi 

and from (3.21) 



= -(2 + 7/3), (3.32), 



^=Z 3 + 1. (3.33) 
From (3.31), (3.33) after excluding L\ and M2 we obtain 

K1M1 = — (2 + jfa)(fa + 1). (3.34) 
As K\ and Mi are positive, and Z3 < — 1, then 

Z 3 > --. (3.35) 

7 

At 7 > 2 it goes already to the contradiction. 

Further, from (3.32) and (3.33) we have K ^l ti — ~ this fact together with (3.34) 
give 

L\ = {2 + 1 faf. (3.36) 
From (3.31), (3.34), (3.35) and (3.36) we obtain the inequality 

-- < -1- v / l-(2 + 7^) 2 -2(2 + 7Z 3 )a3 + l), (3.37) 
7 

which cannot be true at 7 > 1. 
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If h < -1, l 2 = -1, ^3 < -1, then from (3.23) wc get p 2 = -7X4 + 

2L 2 N 



M 2 



= -2-7/3, (3.38) 



from (3.25) q = -7L4 + 



^ = -2-7^3, (3.39) 



from (3.27) 

P2 = --fLi = 7/3. 

In that way, Is = — -, this fact together with (3.38) or (3.39) contradicts to the fact that 
L 2 , N and Mi are not equal to zero. 

If h < -1, h = -1, h = -1, then from (3. 23), (3. 25) and (3.27) we get correspondingly 

2L 2 N L X M X 

+ -71— = -2 - 7^3, (3.40) 



M 2 M 2 
L 2 Mi 



2N 
L X M 2 



from (3.20), (3.21) 



Mi 

K\M 2 2K X N 



= -2-7^3, (3.41) 
= -2-7i 3 , (3.42) 



J 3 + 1. (3.43) 



Li L2 

In that way, multiplying (3.40) by (3.42), taking into account (3.43) we get 
from (3.41), (3.42), (3.43) 

Mi 2 = M|, (3.42) 

that is 

(2 + 7 / 3 ) 2 = L\±L\. (3.45) 

As above, from (3.31), (3.34), (3.35) and (3.45) we get the inequality (3.37) which cannot 
hold at 7 > 1. 

It remains the unique possibility l\ = l 2 = I3 = — 1. In the case p 2 < —2, q < — 2 and, as 
follows from (3.23), p x = p 2 or pi = q, that is p\ < —2. Therefore, as follows from (3.29) 

L\ + L\ - L\ = 1. (3.46) 

If pi = p 2 , then \M 2 \ < Mi. If p\ = p 2 , q < pi, then from (3.23), (3.27) we obtain that 
P2 + 7 = ^j^, Pi + 7 = T7T' hence M i = M 2 ) it goes to the contradiction. 
If pi = q*p 2 < pi, then from (3.23), (3.25), (3.27) we get 
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where A = q + 7 = pi + 7. It follows, in particular, that L 2 = —XL\, L\L 3 = — XL 2 , A 2 = 
—L 2 L 3 . Besides, if G\(t) ~ Nit qi , G2{t) ~ N 2 t q2 , t — > 00, where N\,N 2 are some positive 
constants, gi = g 2 = Pi = <7, ^i = AT 2 , Mi = 2iVi, and N 2 > iV 2 in virtue of G 1 (t)G 2 (t) - 
G\(t) > 0. In such way, from (3.47) we get that L\ > A 2 , A 2 > L 2 , L\ > L\, and taking into 
account (3.46), L\ < 1, L\ < 1, A 2 < 1. That is if A > 1 (7 > 3), the conditions mentioned 
in the paragraph cannot hold together. 

At last, consider the case pi = p 2 = q. Then from (3.23), (3.25), (3.27) we have A = 

^ + T = ^- ^ = T7T-Tff' that is the s y stem of lincar homogeneous 
equations with respect to the variables Mi, M 2 , N 

L1M1 - AM 2 + 2L 2 N = 0, 

L 2 Mi - L 3 M 2 - 2XN = 0, 

AMi - iiM 2 + 2L 3 N = 0. 

For the existence of its nontrivial solution the determinant of the system must be equal to 
zero, i.e. 

\L\ - 2XL 2 L 3 + L\L\ + L\L\ - A 3 = 0. 

Taking into account (3.46), involve the function with respect to the variables L 2 and L 3 , 
where A plays the role of parameter, namely 

vf A (L 2 , L 3 ) = A(l - L 2 + L 2 3 ) - 2XL 2 L 3 + [L\ + L§)0TZ| + L§ - A 3 . 

By the standard methods one can show that the function is not equal to zero at A > 1, i.e. 
at 7 > 3 (remember, that p\ < —2). 

So, it remains to investigate the case 1 < 7 < 3. 

Consider equation (3.11), which can be written as 

A'(i) = 2di(t)A(i). 

If we suppose that A(t) ~ const ■ t m , t — > 00, m is a constant, then m = 2L4 = 2. From 
(3.11) we can get E p (t) ~ const • t 1_7 (^ 0), £ — ► 00. Therefore, if we denote E the quantity 
of the total energy of the system, then Ek(t) = E — E p {t) ~ E(l — Cit 1-7 ), here and further 
Ci are some positive constants. But Ek(t) ~ C 2 Gi(t)t~ 2 A~ L 2~ ~ C 3 GkV~ x , where Gfe is at 
least one of functions Gi, G 2 or G 3 . That is G fe - G 4 (l - Gii 1 - 7 )^- 7 - Ggi 1 " 7 . But then 
1— 7<pi<— 2,7>3, and 7 do not belong to the interval under consideration. 
So the proof of the proposition is over. 

Remark 3.1. One can show more shortly, that in the physical case 1 < 7 < 2 the 
situation l\ = l 2 = l 3 = — 1 if impossible. Taking into account (3.10) we have GiG 2 — G\ = 
A~ 7 ~ const ■ i~ 27 , t — > 00. But the degree of the leading term of the expression G\G 2 — G\ 
is not greater then 2pi, therefore p\ > —7, as pi < —2, then 7 > 2. 

Remark 3.2. Actually ; 3 = -L 4 = -2, G X G 2 - G§ ~ consi • t~ 47 , t -» 00. 

Remark 3.3. In the case \i > 0, I = the asymptotic A(i) ~ jE, t — > 00 holds as well, 
that is according to Corollary 2.1 the corresponding solution with a linear profile of velocity 
is interior. If u > 0, I 7^ it occurs that A(£) ^ A j, but the matrix A under arbitrary 
initial conditions A(0) has not to be of form SE + U 2 (see the denotations in the statement 
of Corollary 2.2), i.e. the Corollary cannot be applied in the last case. 

One can write out the conditions of compatibility of initial data of the density and the 
pressure for the case of solution with arbitrary linear profile of velocity. 
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